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Variability on time scales St < t is observed in many gamma-ray burst afterglows. It is well known that 
there should be no such variability if the afterglow is emitted by external shock, which is produced by 
the interaction of ultrarelativistic ejecta with the ambient interst ellar medium, w ithin the framework of 
simple models. The corresponding constraints were established by lloka et al.l l|2005l ) and in some cases are 
inconsistent with observations. On the other hand, if the motion is not relativistic, then the fast variability 
of the afterglow can be explained much more easily. 

In this connection we discuss various estimates of the time of the transition to subrelativistic motion in 
GRB source. We point out, that this transition should occur on an observed time scale of ~ 10 days. In 
the case of a higher density of the ambient interstellar medium ~ 10 2 -10 4 cm -3 or dense stellar wind with 
M ~ 10 -5 -10 -4 M0/year the transition to a subrelativistic motion can occur on a time scale of ~ 1 day. 
These densities may well be expected in star-forming regions and around massive Wolf-Rayet stars. 
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INTRODUCTION 



According to current models of gamma-ray bursts 
(GRBs), ultrarelativistic motion of the jet pointing to- 
ward the observer takes p l ace in their sources ( s ee re- 
views bv iMeszarosl 120021: IZhang and Meszarosi 120031: 



Fir an, 120051) . Their afterglows in X-ray, optical, and 



radio bands are explained by the radiation that emerges 
at the front of the external shock produced by jet inter- 
action with the ambient interstellar medium surround- 
ing the source. In this case, a variable light curve 
cannot emerge from a homogeneous spherical emitting 
shell, that moves inside a cone with the opening angle 
9 > 7 _1 , where 7 — the Lorentz factor of the jet. For 
the remote observer, the variability time scale should 
be St ~ t. 

However, the variability of optical afterglows at 
time scales 5t < t is observed in many cases. The 
best-known is GRB 030329 whose afterglow was stud- 
i ed in great d e tail due to its ex c eptional brightness 
dBurenin et all 120031: lUrata et all . 12004 iLipkin et al . 
2004). Similar variability was also observed in other 



cases where detailed measurements of the light curve 



were obtained: 



2003; Bersier et al 



GRB 02100 4 
20031 



see e.g. 



Holland et al 



do Ugarte Pos tigo et al 



(|Dai et all 120071 iKhamitov et al.l . l2007h and others 



20051). 050 4 08 (Ide Ugarte Postigo et al.| . l2007t)T 060526 



Various explanations for the vari ability of GRB af - 
terglows were discussed in detail by Ioka et al. ( 20051 ). 
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who also established constraints on the afterglow vari- 
ability time scales and amplitudes for various models. 
They are all based on the assumption of ultrarelativis- 
tic motion of the jet. In some observed afterglows, the 
constra ints obtained in simple models of t he jet are vi- 
olated (jloka et all 120051 IKhamitov et~ail 120071 ) . 

However, the motion may well become moderately 
relativistic on the observed time scales in case of higher 
density of the ambient medi um or in presence of a dense 
stellar wind. For example, Dai fc Lu ( 19991 ) discussed 
this possibility for the afterglow of GRB 990123 and 
showed that the transition to a nonrelativistic motion 
on a time scale of about 2.5 days occurs if the density of 
the surrounding medium is ~10 6 cm -3 . In this paper, 
the sideways exp ansion o f the jet once the condition 7 < 
9- 1 is satisfied (jRhoadsl . Il997t 11999b ISari et all Il99flh 
was not taken into account. Because of this expansion, 
the jet collects more material on its way and slows down 
faster. 

In this note, we discuss various estimates of the time 
of the transition to subrelativistic motion and show that 
it may well be ~ 1 day. 



WHEN DOES THE ULTRARELATIVISTIC 
MOTION END? 

The most simple estimate of the time pf the tran- 
sition to nonrelativistic motion, ijvit, taking into ac- 
count the sideways expansion of the jet, was obtained 
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bv lWaxman et all (| 19981) : 

t NR « 27(^2 M) 1/3 0o 7 i 3 days, 



(1) 



where -E52 is the "isotropic" energy of the shell in units 
of 10 52 erg/s/cm 2 , n\ is the particle number density in 
the medium surrounding the source in units of cm , 
and #o.i = 0/0.1 — is the jet opening angle. 

This estimate was obtained using the dependency 
7(£) for adiabatic shell and assuming that the motion 
after the onset of jet sideways expansion approaches a 
spherically symmetric one and take place according to 
the same adiabatic solutions. In order of magnitude, 
it agrees well with the other estimates that can be ob- 
tained from a more comprehensive analysis of the jet 
dynamics. For example, if the jet begins to expand at 
7 = 9 , then, given that after that the gam ma factor 
depen ds on the observed time as 7 cx V 1 ! 2 ( Rhoadsl . 
1999), we will obtain: 



tNR — tj 



(2) 



where tj is the time of the beginning of jet si deways ex- 
pansio n. If we take the expression for tj from Sari et all 
(Il999l) . ;hen tj^R = 26 days with the same parameter 
d ependency as a bove. 

iRhoadd (|l999h assumed that the jet in its rest frame 
expands with the speed of sound c s — c/s/3, not with 
the speed of light. Accordingly, the jet expansion be- 
gins later, at 7 w (3v20) -1 . Using the j(t) depen- 
dency from this work and applying the same inference 
as abo ve, we obtain tNB. = 25 days for th e same param- 
eters. IPanaitescu and Meszaros (|l999l) assumed that 
the jet also expands with the speed of sound. They 
pointed out that, when 6 > 0.1, the shell ceases to 
be ultrarelativistic even before the onset of the side- 
ways expansion. In this work, tNR was estimated ex- 
plicitly. Rewriting it with our parameters, we obtain 
tNR, — 6.1 days, with the same parameter dependency 
as in ([TJ) . 

If we take into account that initially the motion of the 
shell may not be adiabatic and it can lose a substantial 
part of its energy through radiation, then the estimate 
of t n r can decrease significantly. For synchrotron ra- 
diation, the motion should be radiative if much of the 
energy that the material gains at the shock front goes 
into the accelerated electrons and if the electrons cool 
down rapidly compared to the dynamical time. The 
latter always holds at the beginning of shell motion. 
The expression for the energy, to be used to calculate 
the subsequent adiabatic m otion if the evolution was 
initially radiative is given bv lSari et al.l (|1998l ). Substi- 
tuting this energy into ((T|) yields: 



tNR,; 



7.6 e, 



1/5-1/5 p4/15 -4/15 -7/15.2/3 



3 £ 5 T 72 



0Q.i da y s > 



where €r and e e are the energy density fractions of 
the magnetic field and accelerated electrons behind the 
shock, 72 = 70/IOO is the initial gamma factor of the 
ejecta. This estimate should be used if it is assumed 
that e e ~ 1, i.e., the electrons are accelerated effectively 
at the shock front. 

The estimate of t^R for a stellar wind with 
a density proportional t o r~ 2 can be taken from 
Livio and Waxman (120001) . Rewriting it with our typi- 



cal parameters, we will obtain: 

t NR , w « 5.7 E^{M^/v s )- l 6l s days, 



(4) 



where M—5 is the mass loss rate of the star in units 
of 10~ 5 MQyear" 1 and W3 is the wind velocity in units 
10 3 km s _1 . If the electrons are accelerated effectively 
and the evolution is initially radiative, then this esti- 
mate should be much smaller, just as in the above case 
of a constant-density medium. 

DISCUSSION 

From these estimates we see, that the observed tran- 
sition time to nonrclativistic motion is ~ 10 days even 
for the commonly accepted typical parameters. A 
higher density of the ambient medium n ~ 10 4 cm -3 for 
an adiabatic shock, or n ~ 10 2 cm -3 , if the electrons at 
the shock front are accelerated effectively and the shell 
is initially radiative, or a stellar wind with M ~ 10 -5 - 
10~ 4 Mq year -1 is required for t^R ~ 1 day. 

A higher density of the interstellar medium is ac- 
tually may well be expected near GRB sources, since 
at least a large fraction of them are related to super- 
nova explosions at the end of the evolution of massive 
stars and occur in regi ons of enhanced star forma tion 
(see, e.g., the review bv lWooslev and Bloo m. 20061 and 
references therein). In addition, massive stars at the 
end of their evolution intensively lose their outer lay- 
ers a nd should be su rrounded by dense stellar winds 
fe.g.. lCrowthen . l2007l ). In some cases, the presence of 
a wind around GRB sources is confirmed observation- 
ally For example, high-resolution spectroscopy of the 
GRB 021004 afterglow shows that its progenitor was a 
massive Wolf-Rayet star surrounded by a wind with a 
mass loss rate of ~ 10~ 4 Mq year " 1 and an the e xpan- 



sion velocity up to 3000 km s 1 ( Mirabal et al. . 2003t 



Lazzati et all . l2006h 



(3) 



If the motion is no longer ultrarelativistic, then the 
fast variability of the afterglow on a time scale St <C t 
can be explained, for example, by the presence of den- 
sity inhomogeneities. In fact, the presence of these in- 
homogeneities can also be expected. For example, the 
stellar winds around Wolf-Rayet stars are known to be 
highly inhomogeneous and to have a clumpy s tructure 
([Hamann and Koesterkall998uCrowtheill2007l ). In ad- 
dition, the cooling time in the observed optical part of 
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the spectrum should be short. This holds for the syn- 
chrotron sp ectrum at late e volutionary phases of the 



shock (e.g. JSari et all 119981 ). 



Thus, the transition to subrelativistic motion on the 
observed time scale of ~ 1 day is actually quite possible 
and, in this case, the fast variability of GRB afterglows 
can be easily explained. Of course, the fast variability 
of GRB afterglows can also be explained under spe- 
cial assumptions in the case of an ultrarelativistic jet 
( Ioka et al l l2005l) . however, the transition to a nonrel- 



ativistic motion should at least be considered as one of 
the possible explanations of this fast variability. This 
requires a higher density interstellar medium or dense 
stellar wind, which are expected in star forming regions 
and around massive Wolf-Rayet stars. 

The author is grateful to the anonymous referee 
who made several important remarks regarding the 
content of the paper. This work was supported by 
the grants RFFI 05-02-16540, RFFI 07-02-01004, NSh- 
1100.2006.2, MK- 4064.2005.2, and Programs P-04 and 
OFN-17 of the Russian Academy of Sciences. 



P. Meszaros, Ann. Rev. of Astron. and Astrophys. 40, 
137 (2002). 

Mirabal N., Halpern J. P.,Chornock R., et al, Astro- 
phys. J. 595, 935 (2003). 

Panaitescu A., Meszaros P., Astrophys. J. 526, 707 
(1999). 

Piran T., Rev. Mod. Phys. 76, 1143 (2005). 

Rhoads J. E., Astrophys. J. (Letters) 487, LI (1997). 

Rhoads J. E., Astrophys. J. 525, 737 (1999). 

Sari R., Piran T., Narayan R., Astrophys. J. (Letters) 

497, L17 (1998). 
Sari R., Piran T., Halpern J. P., Astrophys. J. (Letters) 

519, L17 (1999). 
Urata Y.,Miyata T., Nishiura S., et al, Astrophys. J. 

(Letters) 601, L17 (2004). 
Waxman E., Kulkarni S. R., Frail D. A., Astrophys. J. 

497, 288 (1998). 
Woosley S. E., Bloom J. S., Ann. Rev. of Astron. and 

Astrophys. 44, 507 (2006). 
Zhang B., P. Meszaros, International Journal of Modern 

Physics A 19, 2385 (2003). 



REFERENCES 
Bersier D., Stanek K. Z., Winn J. N., et al., Astrophys. 

J. (Letters) 584, L43 (2003). 
Burenin R. A., Sunyaev R. A., Pavlinsky M. N., et al., 

Astronomy Letters 29, 573 (2003) [Pis'ma v Astron. 

Zhurn., 29, 649; |astro-ph/0306137| . 
Crowther P. A., Ann. Rev. of Astron. and Astrophys., 

in press (2007); |astro-ph/0610356| 
Dai Z. G., Lu T., Astrophys. J. (Letters) 519, L155 

(1999). 

Dai X., Halpern J. P., Morgan N. D., et al., Astrophys. 

J. 658, 509 (2007). 
de Ugarte Postigo A., Castro-Tirado A. J., Goros- 

abel J., et al, Astron. Astrophys. 443, 841 (2005). 
de Ugarte Postigo A., Fatkhullin T. A., Johannesson G., 

et al, Astron. Astrophys. (Letters) 462, L57 (2007). 
Hamann W.-R., Koesterke L., Astron. Astrophys. 335, 

1003 (1998). 

Holland S. T., Weidinger M., Fynbo J. P. U., et al, 

Astron. J. 125, 2291 (2003). 
Ioka K., Kobayashi S., Zhang B., Astrophys. J. 631, 

429 (2005). 

Khamitov I. M., Burenin R. A., Bikmaev I. F., et al, 
Astronomy Letters 33, 797 (2007) [Pis'ma v Astron. 
Zhurn., 33, 891]. 

Lazzati D., Perna R., Flasher J., et al, Mon. Not. R. 
Astron. Soc. 372, 1791 (2006). 

Livio M., Waxman E., Astrophys. J. 538, 187 (2000). 

Lipkin Y. M.,Ofek E. O., Gal- Yam A., et al, Astro- 
phys. J. 606, 381 (2004). 



ASTRONOMY LETTERS Vol. 33 No. 12 2007 



